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I ^ ■ Abstract 

I , We present the full set of analytical solutions of the geodesic equations of charged test particles in 

the Reissner-Nordstrom space-time in terms of the Weierstrafi p, a and £ elliptic functions. Based on 
the study of the polynomials in the i? and r equations we characterize the motion of test particles and 
discuss their properties. The motion of charged test particles in the Reissner-Nordstrom space-time 
' ' ■ is compared with the motion of neutral test particles in the field of a gravitomagnetic monopole. 

^ji Electrically or magnetically charged particles in the Reissner-Nordstrom space-time with magnetic 

or electric charges, respectively, move on cones similar to neutral test particles in the Taub-NUT 
^q' space-times. 



The Reissner-Nordstrom space-time is a static, asymptotically flat solution of the Einstein-Maxwell equa- 



1 Introduction 

> 

tions in general relativity. It describes charged, non-rotating spherical black holes or naked singularities. 
■ In general, the Reissner-Nordstrom space-time represents a gravitating source which is both electrically 

and magnetically charged. It is a special case of the general family of electrovacuum space-times of Petrov 
type D found by Plebahski and Demiahski [T] and reviewed by Griffiths and Podolsky [2J. These space- 
times describe the gravitational fields of isolated massive objects, e.g. stars or black holes, and are in 
addition characterized by a NUT charge, angular momentum, acceleration and a cosmological constant. 
The higher dimensional generalizations of the Reissner-Nordstrom solution were found by Tangherlini [3] . 

The physical properties of a space-time can be investigated by studying the motion of test particles 
and light in this space-time. Test particles are a sensitive tool to investigate the properties of the fields 
generated by the massive body because of the coupling of the parameters of the metric and the test particle 
in the equations of motion. The geodesic equations can be solved either analytically or numerically. A 
pioneer in finding analytical solutions was Hagihara [4]. He integrated the geodesic equation of a test 
particle in the Schwarzschild gravitational field in terms of the Weierstrass p-function. Some special cases 
of the geodesies in the Reissner-Nordstrom space-time were studied in [SJ [3]. General features of the 
motion of charged test particles in the Kerr- Newman space-time were considered in [7J. The influence 
of the gravitomagnetic monopole moment on the motion of test particles in Kerr-Newman- Taub-NUT 
space-time was studied in [5]. 

A breakthrough in the analytical integration of geodesies came with the papers of Hackmann and 
Lammerzahl [5J [TU] where geodesies in a 4-dimensional Schwarzschild-de Sitter space-time were integrated 
analytically in terms of the hyperelliptic 9 and a functions. The developed method is based on the 
Jacobi inversion problem restricted to the ^-divisor. Also geodesies in the Reissner-Nordstrom-de Sitter 
space-time in 4 dimensions as well as in some higher dimensional Schwarzschild, Schwarzschild-(anti)de 
Sitter, Reissner-Nordstrom and Reissner-Nordstr6m-(anti)de Sitter space-times were integrated by this 
mathematical method |llj . The elliptic and hyperelliptic functions were subsequently used to obtain the 
solutions of the geodesic equations in the axially symmetric Taub-NUT (T2] and Kerr-de Sitter [T3] space- 
times. In these papers the types of orbits are classified and the orbits of test particles are extensively 
studied. Furthermore, the analytical solution of the geodesic equations in the general Plebahski-Demiahski 
space-time in 4 dimensions were obtained in terms of the hyperelliptic Kleinian a function |14) . 

The Reissner-Nordstrom space-time as well as the Schwarzschild space-time have singularities at their 
origin. Whereas in the Schwarzschild space-time particles are bound to end up in the central singularity, 



this does not hold for particles which cross the event horizon in the charged space-time. Here the potential 
barrier due to the charge does not allow the singularity to swallow the particles. Instead, a (test) particle 
will leave the vicinity of the singularity again. After traversing the Cauchy horizon and the event horizon 
of the Reissner- Nordstrom black hole, the particle will emerge in another universe [5j [15l [IT]. Moreover, 
an orbit in such a black hole space-time can be analytically continued into an infinite sequence of further 
patches of the space-time. This feature of the Reissner-Nordstrom solution is seen clearly in its Carter- 
Penrose diagram. 

The Taub-NUT space-time which is characterized by a gravitomagnetic charge possesses also two 
horizons. A Kruskal-like analytic extension [16] of the Taub-NUT space-time would then likewise allow 
test particles to move into other worlds. However, such an extension comes at the prize that the periodic 
identification of the time coordinate, suggested by Misner to remove the singularity from the symmetry 
axis, is no longer possible [16j . Consequently, the space-time is geodesically incomplete, since orbits end 
either at the singular axis (in the Bonnor-Manko-Ruiz interpretation) or at one of the horizons (in the 
Misner- Taub interpretation). For a recent discussion see [T2"] . 

But the presence of two horizons and the associated many-world orbits is not the only similarity 
between the Reissner-Nordstrom space-time with electromagnetic charge and the Taub-NUT space-time 
with gravitomagnetic charge. Also the motion of test particles has common properties in both types of 
space-times. In this paper we study the motion of electrically and/or magnetically charged particles in 
the general Reissner-Nordstrom space-time. As in the Taub-NUT spacetime where a test particle moves 
on a cone, the motion of a charged test particle in the general Reissner-Nordstrom space-time proceeds 
on a cone. The condition for the motion on a cone is that the products of the dual charges should not be 
equal i.e., Gq 7^ Qg (where Q and G are the electric and magnetic charges of the gravitational source and 
q and g those of the test particle, see Section [5] for details). Otherwise the motion would be planar. 

In the following section we present the equations ot motion. We classify the orbits of charged test 
particles in Section [3] We analyze their motion by studying the influence of the parameters of the metric 
and of the test particles. In particular, we consider the dependence on the mass, on the angular momentum, 
on the charges and on the separation (Carter) constant. In Section 0] we solve analytically the r, i9, (f 
and t equations of motion, and we give the solutions of the radial and time equations in terms of the 
Weierstrass elliptic p, a and £ functions. We give details on the elliptic integrals in the Appendices. 



2 The geodesic equation 



The Reissner-Nordstrom solution of the Einstein-Maxwell field equations is described by the metric [2] 

ds 2 = ^-dt 2 - ^-dr 2 - r 2 idd 2 + sin 2 tfdtp 2 ) , (1) 

where A r = r 2 — 2Mr + Q 2 + G 2 . Here M is the mass of the solution, Q and G are the electric and 
magnetic charges. The singularity is located at r = 0. For < Q 2 + G 2 < M 2 , there are two horizons, 
defined by A r =0, and given by 



r ± =M(l±y/M 2 -{Q 2 + G 2 )) . (2) 



Between the horizons the radial coordinate r becomes timelike, and the time coordinate t spacelike. When 
Q 2 + G 2 = Al 2 , the horizons degenerate and the black hole is called extremal. For Q 2 + G 2 > M 2 the 
solution has a naked singularity. 

The field strength = A v ^ — A^^ and the dual field strength F^ IV = A V4l — A^^ of the electromag- 
netic field are induced by the non- vanishing components of the vector potentials A^ and A^ 

A t = -, A tp = -Gcosti, (3) 
G 

A t =i— , ^L = iQcosz? . (4) 
r 

The dual field strength is defined by the antisymmetric Levi-Civita symbol e^ vaT as F^ v = — * £ a«"" ' F aT 

\ 9 

with g d = -det||g M1/ ||. 



The Hamilton- Jacobi equation for a particle with electric charge q and magnetic charge g 



dS uv ( dS j \ / OS , \ 

2 d^ = 9 \d& qA » + igA » ) [d^ qK + igAv ) (5) 

has a solution in the form S — \5t — Et + Lip + S(r) + S($). Here r is an affine parameter along the 
geodesic, S is a parameter which is equal to for a massless test particle and equal to 1 for a test particle 
with non-zero mass. The constants E and L are the conserved energy and the angular momentum in the 
z direction of a test particle, respectively. Because of the presence of the angle dependent terms qA v and 
gA v in ((5|) the problem of charged particle motion is axially symmetric, although the space-time itself as 
described by the metric ([1]) is spherically symmetric. 

For convenience, we introduce dimensionless quantities (rs = 2M) 

~ r ~ t ^ t ~ Q ~ G ~ L 

r = — , t= — , t = — , Q = — , G = — , L = — . (6) 

rs r s r s r s r s r s 

The Hamilton- Jacobi equation ((5} separates and yields for each coordinate a corresponding differential 
equation 

df s 2 

dd 



- R (7) 
= e (8) 

dip 1 



L + A g costf (9) 



dj sin 2 1? 



with the polynomial R and the function O 



R = f i [E+^\ -K r {5r 2 + k) (11) 



9 = k 2— L + Agcostf . (12) 

sin d V / 

The following notations were introduced: A g = Gq — Qg, A q = Qq + Gg and A r = f 2 — f + Q 2 + G 2 . Here 
k = K + L 2 , where K is known as Carter constant or separation constant. For zero charges the equations 
of motion above reduce to the Schwarzschild case. We also used the Mino time 7 as r 2 d^ — df [17]. Note, 
that for A g = the situation simplifies, and the motion takes place in a plane. 



3 Complete classification of geodesies 

Consider the Hamilton- Jacobi equations (|7|)- (fT0|) . The properties of the orbits are given by the polynomial 
R (|11[) and the function (|12[) . The constants of motion (energy, angular momentum and separation 
constant) as well as the parameters of the metric and the charges of the test particle characterize these 
polynomials and, as a consequence, the types of orbits. In this section we discuss the charged motion in 
the Reissner-Nordstrom space-times in terms of the properties of the underlying polynomial R and the 
function 6. 



3.1 The motion 

Since $ is a polar angle, it can take only real values. Equation © has real solutions if the requirement 
> is fulfilled. This implies k > 0. With the new variable £ = cos Eq. © turns into the equation 



dC 2 



= 65 with 0£ = < 2 + + c , (13) 



id o 

c = k — L 2 . From k > follows a < 0. The zeros of define the angles of two cones which confine 



with a simple polynomial of second order on the right hand side, where a = —(k + A 2 ), b = — 2LA g , and 

"1 



the motion of the test particles (a similar feature appears in Taub-NUT [T2J and Kerr space-times [T5]). 
Moreover, every trajectory is not only constrained by these cones but lies itself on a cone in 3-space [T81IT5] . 
If A g in the equation Q and in the polynomial (|12p vanishes, then the motion lies on a plane (e.g., the 
motion of only electrically charged or neutral particles reduces to a plane in a Reissner-Nordstrom space- 
time with only electric charge). In the space-time of a gravitomagnetic monopole the trajectories of test 
particles similarly lie on cones [201 [5TJ [55] . 

The discriminant D = b 2 — Aac of the polynomial can be written as D = 4kn with k = k + A 2 — L 2 . 
The existence of real zeros of 0£ requires D > 0. This implies that both k and k should be non-negative 



k > 

K = k- L 2 + A 2 > 



2 (14) 



These are conditions on the parameters L and k for given values of Q, G, q and g. As long as k — L 2 is 
positive there are no constraints on L. When k — L 2 becomes negative the inequalities (|14p imply a lower 

limit for the angular momentum given by L m i n = ±y/A 2 + k. 

One can show that 



lf~ \ 2 1 /- \ 2 

© = fc g— (Z + costfAg =k „— (Zcostf + Ag . (15) 

sin ■& V / sin ft V / 

The zeros of ©£ are given by 



_ LA g ± Vkn 

a ' 2 Hk + Aj) ■ (16) 

The conditions (fl4| ensure the compatibility of £ G [— 1, 1] with > 0. 

/ LA„ kn \ 

The function Qc describes a parabola with the maximum at — — , — . For non- vanishing 

4 y k + Aj' k + AjJ 6 

L and A g the maximum of the parabola is no longer located at £ = or, equivalently, the zeros are no 
longer symmetric with respect to £ = 0. Only for vanishing L or A g both cones are symmetric with 
respect to the equatorial plane. 

The i9-motion can be classified according to the sign of k — L 2 : 

1. If k — L 2 < then Oj has 2 positive zeros for LA g < and $ £ (0,7r/2), so that the particle moves 
above the equatorial plane without crossing it. If LA g > then d 6 (ir/2,ir). 

2. If k - L 2 = then 6< has two zeros: & = and £, = - - 2 2 £A A 8 , . If LA g < then £ e [0, 1) and 

i9 € (0, §]. If LA g > then £ G (-1,0] and G [f ,tt). If Z = -A g then the ^-motion fills the 
whole upper hemisphere i? G [0,5]. The motion fills the whole lower hemisphere with ■d G [§,7r] if 
X = A g . 

3. If k — L 2 > then has a positive and a negative zero and •& G (0, 71"), and the particle crosses the 
equatorial plane during its motion. 

In general, the second term of the function O in (fT5|) diverges for d — > 0, ir. However, if L — — A g this 
term is regular for $ — and if L = A g it is regular for i? = tt. If L = ±A g , then k — n. The regularity 
of O in these cases can be seen from 



?2 (lTco_s_tf) 
sin J 



2 



e = K -L ^ \,__ 2a > , (17) 



by application of L'Hopital's rule. If, furthermore, L 2 = then O = k + Ag — L 2 which is independent of 
<&. 

In the special cases when one of the constants k or k or both vanish, O^ has a double root which is 
the only possible value for §. One can distinguish three cases: 

1. If k = and k > 0, then £ = for L 2 < A 2 

2. If k = and k > then £ = ~t* for i 2 > A 2 



3. If k = k = then £ = ±1 implying that d — or i3 = it are possible. In this case L = T^g (as 
discussed above). 

This means that during a test particle's motion the coordinate i9 is constant and the trajectory lies on a 
cone around the i9 = 0,7r-axis with the opening angle arccos£. In this case we immediately obtain from 
© that 93(7) = C(-/ — 7in) with a constant C = -jrgs 3 which in case 1 is C = and in case 2 is C = L. 

Thus, the non-vanishing of k and K indicates that the motion of the particle is not symmetric with 
respect to the § = 0, 7r-axis. Therefore, these two constants may be regarded to play the role of a 
generalized Carter constant which appears, e.g., in the motion of particles in NUT and Kerr space-times. 

3.2 The f— motion 

3.2.1 Possible types of orbits 

In the following we consider the motion of charged particles in the regular Reissner-Nordstrdm black 
hole space-time, possessing two non-degenerate horizons r±. Before discussing the r-motion in detail, we 
introduce a list of all possible orbits: 

1. Escape orbits (EO) with ranges: (n, 00) with r% > r + . These escape orbits do not cross the horizons. 

2. Two-world escape orbits (TEO) with ranges (ri,oo) with T\ < r^ 

3. Periodic bound orbits (BO) with range f € (r%,r2) with T\ < r2 and 

(a) either n , r 2 > r + or 

(b) n,r2 < r_. 

4. Many-world periodic bound orbits (MBO) with range f £ (r\,r2) where r% < r_ and r2 > r + . 
Here we used the notation for the orbits introduced in [11] for the regular Reissner-Nordstrom space-time. 

3.2.2 The radial motion 

The right hand side of the differential equation ([7]) has the form R — Ylt=o with the coefficients 

64 = E 2 -5 (18) 

b 3 = 5 + 2EA q (19) 

b 2 = -(k- Al + S{Q 2 + G 2 )) (20) 

61 = k (21) 

b = -k(Q 2 +G 2 ). (22) 

Let us now consider massive particles only, that is 8 = 1. In order to obtain real values for f from (J7J we 
have to require R > 0. The regions for which R > are bounded by the zeros of R. The number of zeros 
depends on the values of E, k, Q, G, q and g. Two conditions R = and 45 = define the double zeros 
of the polynomial R and thereby the boundary between the regions where R has 1, 2, 3 or 4 zeros. The 
parameter plots shown in Fig. [1] are based on this. One has to additionally take care of the change of the 
sign of E 2 — 1 when E crosses E 2 = 1. Then the sign of R(r) for f —> ±00 changes. Furthermore, the case 
E 2 = 1 requires additional attention: If the line E 2 = 1 is contained in a region with 3 or 1 zeros then on 
this portion of the line we have 3 or 1 zeros respectively. Taking all these features into account we obtain 
the k-E diagrams of Fig. Q] 

We define the effective potential from equation (JT)) as the values of energy E when 

df\ 2 

=r*(E-V+)(E-V cS ), (23) 
thus, 

y c ± S = -^ L± ^\/^r(Sf 2 + k). (24) 

Condition (|23l) then determines the turning points of an orbit. 

Examples for the effective potential are given in Figs. [2l [3l These figures combine the positive-root 
part V^g and the negative-root part V^. At the horizons both parts coincide (i.e. at these points 



both parts are glued together), since A r vanishes here. Consequently, V^ f (r±) = — Clearly, 

Kff( r ±) < if A q > and vice versa (see e.g. Fig|3|). 

We note, that in contrast to the Schwarzschild case here the positive-root V^g can allow for particles 
with negative energies. The region outside the horizon where negative energy particles may sojourn has 
been termed 'generalized ergospere', since energy may be extracted [23j [24]. While particle energies 
below the negative-root have no classical interpretation, they are associated with antiparticles in the 
framework of quantum field theory [25]. The relation V^(E,q,g) = —V^{—E,—q,—g) shows, that the 
corresponding positive energy orbits are available for particles with opposite charges (i.e., antiparticles). 
A lower limit for the energies of particles is obtained from the requirement that time should only run 
forward (see Eq. ([TO])). E > -fe. 

As we know from previous studies for neutral particles (see e.g. [TT]), the Reissner-Norstrom space- 
time possesses a potential barrier which prevents particles from falling into the singularity. The potential 
barrier which is defined by the smallest positive root n of R is located in the interval < n < r_ . 

For increasing values of A q , the influence of the charge of a test particle becomes noticeable through 
the term (We recall, that the charges of the black hole Q and G are constrained by the naked 

singularity condition Q 2 + G 2 > i). Namely, for larger values of the test particle charge the effective 
potential can form a small potential mound (Fig 3(a) I or potential well (Fig |3(b)] ) in the interval (0, \ 



, r_ 



This indicates that bound orbits may exist in this region. Such a feature is only present in the motion of 
charged particles. 

Using the k-E diagrams in Fig. [1] as well as the above considerations we can give all possible combi- 
nations of zeros of R and an interpretation in terms of specific types of orbits which are summarized in 
Tabled! 

The types of orbits related to the various parameters are given by: 

• Region (1): one positive zero. The orbit is a TEO with particles coming from r = +oo. From the 
features of the effective potential (potential barrier) it is clear that the turning point can lie on the 
inner horizon (case A-). In this case the energy of a test particle corresponds to the glue point of 
the potentials V^: Ea_ = (?"*_). For E 2 = 1 the coefficient of the highest power in R, given by 
E 2 — 1 for massive test particles, vanishes. On the line E 2 = 1 within the region (1) there is one 
positive zero. 

• Region (2): two positive zeros. Here only MBOs are possible. In the special case when A q = and 
E 2 = 0, the turning points are lying on the horizons (case B±). Also the cases B- and B+ when 
one of the turning points lies on the horizons are possible. In this case the energy again corresponds 
to the glue points of the potentials V^: Eb_ — V^(r_) (such an orbit can be found in e.g. Fig 2(a) I 
or Eb + = Kff( r +) (see e.g. Fig |3(b)[ ) respectively. 

• Region (3): three positive zeros. 

— Region (3)+: Here MBOs and EOs with corresponding subcases C_ and C+ are possible. 

— Region (3)_: For growing charges the term — -=r leads to a bound orbit with turning points 
behind the inner horizon. It is also possible that a turning point of the TEO lies on the inner 
horizon (case D-). 

Here too for E 2 = 1 the term of highest power in R vanishes giving three positive zeros for the 
E 2 — 1-line within the region (3). 

• Region (4): four positive zeros. Here we find MBOs and BOs (planetary orbits). 



4 Solution of the geodesic equation 

Now we present the analytical solutions of the differential equations (f7l)- (fl"0|) . 
4.1 Solution of the equation 

The solution of Eq. (IT3)) with a < and D > is given by the elementary function 

0( 7 ) = arccos(i (Vd sin ( - it) ~ &)) » (25) 



type region +zeros 



range of f 



orbit 



A (1) 1 • 1 i TEO 

A_ I II TEO_ 

B (2) 2 * II II ■ MBO 

B+ H ♦ * MBO+ 

B_ ♦ b* * MBO_ 

B± I ♦ ' MBO± 

C (3)^ 3 * II II * * MBO, EO 

C_ ♦ II * * * MBO, EO 

C+ *H f • MBO+, EO 

D (3)_ 3 •"H 1 BO, TEO 

D_ * * f ! ' BO, TEO 

E (4) 4 * II II * * * : MBO, BO 



Tabic 1: Types of polynomials and orbits of charged particles in the Rcissncr-Nordstrom space-time. The thick lines represent the 
range of the orbits. The turning points arc shown by thick dots. The horizons arc indicated by a vertical double line. In 
special cases the turning points lie on the horizons. Type D with a bound orbit behind the inner horizon is possible for 
charged particles with relatively large charges. 



where 7 m = a^- ln — arcsin ( ^1'"^, ) an< i 7in is the initial value of 7. 



4.2 Solution of the f— equation 

For timelike geodesies the polynomial R in (fTTj) is of fourth order. A standard substitution f = ±A + fn, 
where fn is a zero of i?, reduces J7| to a differential equation with a third order polynomial [jj^j — 
i?3, where R3 — Y^i=o^ iX% - further substitution x = (4y — ^.) transforms that into the standard 



Weierstrafi form 



d'y 

where 



4y 3 -522/-.93 =Ps(y), (26) 



= 6| _ 6163 _ 616263 _ &o&f _ _6J_ 

92 12 4 ' 33 48 16 216 ' 1 ' 

The differential equation ([25)1 is of elliptic type and is solved by the Weierstrafi p-function [35] 

2/(7) = P (7- 7in5 52,ff3) , (28) 

where 7/ = 7i n + f°° 1 dy with y in = ±^ (r in — f^) -1 + l§. Then the solution of acquires 
the form 

f = ± ^ — r 1 — r^ + ?R - (29) 

4 P (7-7 in ;52,ff3) - f 
4.3 Solution of the (/^-equation 

Eq. ([9]) can be simplified by using (jHJ and by performing the substitution £ = cos 1? 

^ ^ 1 - e 2 ^ 1 - e 

where ©£ is given in (fT3"]) , This equation can be easily integrated and the solution for a < and D > is 
given by 

^(7) = ^(/ + + /-)£ (7) +¥>m, (31) 




(c) Q = 0.3, G = 0.1, q = 10, g = (d) Q = 0.4, G = 0.2, q = -3,3 = 0.1 



Figure 1: Parametric fc-S digrams showing the location of the regions (1)— (4), which reflect the number of zeros of the polynomial R 
in Eq. JlH ■ Each region contains a set of orbits peculiar to it which are described in Tableland in the text (Section [3~X2l . 
Here the influence of the variation of the el ectric ch arge q of a test particle on the zeros of R is presented. Increasing the 
positive value of the electric charge in plots [(a)|(c)| one observes, that the region (4) with four positive zeros for negative 
E (left side) disappears slowly. When th e sig n of q is changed (along with the sign of g) the right and left side of the 
k-E diagrams are mirrored (compare plot |(d)| with the rest). One finds similar k-E diagrams by varying the value of the 
magnetic charge g. Because of the inequalities ]14t the constant k — K + L is positive. 



where 



L±A k + K-(L±A g y=f[k + K + (L±A fs y)C 

I± = — , arcsin =^ '— (32) 

|L±A g | v^(CTl) 

For the special case k = L 2 and L = ±A g the solution reduces to the simple form 



, , 1 L . 1±3£ 
<P\"t) = — =~ arcsm 



2|L| CTl 



+ Vin ■ (33) 



4. 4 Solution of the ^-equation 

Eq. ([10]) consists of two f-integrals: 



Consider I}. The substitution r — ± — + rjj reexpresses If in terms of y: 

4 v-^ 



dy (r R (4y- b f)±b3 4 



/1( 7 ) = E I ^VM^ i^I£7 (35) 
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(a) Q = 0.4 , G = 0.25 , q = 0.05 , g = 0.1 
Effective potential with four cxtrcma. 




(b) Detail of figure (a): 1 to 4 turning points arc possible. 

for orbits of type A, B, C, E from Table [T1 showing 
the position of the orbit types BO, MBO, EO, TEO. 
The points indicate the turning points of the motion 
and the red horizontal dashed lines correspond to the 
values of the energy E. 




(c) Q = 0.4 , G = 0.2 , q = -3 , g = 0.1 , k = 1. 



(d) Q = 0.4 , G = 0.2 , q = 3 , g = 0.1 , k = 1. 



Figure 2: Effective potential Vli(f) for different values of Q, G, q, g and k. The blue line represents V1~L , the green line represents 



V cff . The two potential parts glue at the horizons— glue points— with the ordinate V^(r±) — 



The grey area marks 



a physically forbidden zone. The positions of the horizons arc shown by vertical dashed lines. At infinity the effective 



potential tends to the limiting values linif 



Figures | (c) | and | (d) | arc plotted for opposite electric charges of 



the test particle. The gl ue po ints of the V^ ff arc cither positive (for q < 0) or nega tive (for q > 0). As can be seen from 
the k-E diagram in Fig. |l(d)| there are regions with 1, 2 or 3 zeros for the potential |(c)| with k — 1. 




1,5 




a) Q = 0.4 , G = 0.25 , q = —4 , g — 0.1 , k — 0.2 



(b) Q = 0.4 , G = 0.25 , q = 4 , g = 0.1 , k = 0.2 



Figure 3: Effective potential V cff (f) for different values of Q, G, q, g and k. The grey area marks a physically forbidden zone. In 
this cas e up to 3 turning points exist for f < f _ . Such an effect can be caused by a high charge of the test parcticle. In 
the plot |(a)| the orbits of type A, C, D from Tableware shown. 



where A y = A r (f R ) (Ay - f ) ± (2f R - 1) b 3 (Ay - ^) + &§ = l6A r (r R )(y - Pl )(y-p 2 ). Here, Pl and p 2 
are two zeros of A y . 

We next apply a partial fractions decomposition upon Eq. (|35|) 



4„ 71) 



x 4 



(36) 



where £>3 = j§ anc ^ * = 0^ ■-; 4, are constants which arise from the partial fractions decomposition. 
These depend on the parameters of the metric and the test particle and on tr. After the substitution 
y — p(v) with p'(v) = (— l) p ^/4p 3 (u) — g 2 p(v) — 33, where p is either or 1 depending on the sign of 
p'(v) in the considered interval and on the branch of the square root, Eq. (|3"o) simplifies to 



)-Pj (p(v)-pa) 



dv . 



(37) 



Here v = v(j) = 7 - 7^ and v ln = v(y ln ). 

The second integral I~ in (134)) can be reduced in a similar way to the necessary form, yielding: 



\dv 



(38) 



where Hi, i = 0, .., 3, are constants which arise from the partial fractions decomposition. 
The final solution takes the form (details can be found in the appendices [A] and [B]) 



a(v - vji) 



r(vi n - Vji) 



K i ( H ^ >/ n , P ( v 3i) , Cr(w-W3i) 

'— -2 C « - W3i) ~ C Win - 1<3i + — V l0 § "7 \ 

OP («3i)) V P ( U 3i ) cr (v in - W 3i ) 



T (-l) p AjiJ ( t ;-« in )+E 



£;7 



H, 



C(Vji)(v - Win) + log 



u(vin - Vji) 



(39) 



where Vji are the poles of the functions (p(v) — Pj) 1 and (p(v) — ps) 2 in (|36| such that p(vji) = pj = 
p(v j2 ) since p(u) is an even elliptic function of order two which assumes every value in the fundamental 



parallelogram with multiplicity two. Here is the Weierstrafi zeta function and <r(v) is the Weierstrafi 

ft +• loci A / P( v 3i) . p"(v 3i )C(v 3i )\ . 

Sigma function \2b ; Ao = — > ^ H = — is a constant. 

4.5 The orbits 

With these analytical results we have found the complete set of orbits for massive charged particles. The 
following figures show solutions of the geodesic equations for the spatial coordinates. We start with small 
charges of a test particle in Figs. 21 [5l [6] The orbits here lie on cones with large opening angle A$ —> n. 
Depending on the energy of the test particle one gets MBO, BO, TEO or EO. In all these figures the 
energy is chosen in the vicinity of the maximum of the effective potential. 

In the Figs. [5J El we increased the charge of the test particle. The orbital cone is most evident 



here. The orbit 7(a) with its x-z projection |(b)| is of type D from Table [T] The test particle orbits the 
singularity at f — 0. Such kinds of orbits are pecular to particles with large charge (and small constant 
k). An example of the effective potential is presented in Fig. [3] As compared to super test particles with 
small mass and small charge, a particle with small mass and large charge is not super any more, but it 
may still be regarded as a test body or test particle [57] . 




Figure 4: MBO and BO with parameters Q = 0.4, G = 0.2, q = 0.1, g = 0.1, k = 3, E = 0.9548, L = 1.0. The sphere in the 3d 
plots shows the horizon f+. Orbits lie on cones with large opening angle A$ — > 7T. 




Figure 5: TEO with parameters Q = 0.4, G = 0.2, q = 1.0, g = 0.1, k = 6.5, E = 1.0039, L = 1.0. The particle crosses both 
horizons, is reflected at the potential barrier and emerges into a second universe, where it moves out to infinity. The orbit 
lies on a cone with large opening angle Ai9 — ^ 7r. 



l - 




Figure 6: MBO and EO with par amete rs Q = 0.4, G = 0.2, q = 1.0, g = 0.1, k = 6.5, E = 1.0038, L = 2.5. An increase of L 
(compare e.g. with Fig. |5(b)|| leads to a decrease of the opening angles of the two cones, which confine the orbit. The 
orbits lie on cones with large opening angle A$ — > 7r 




Figure 7: BO and TEO with parameters Q = 0.4, G = 0.25, q = -4, g = 0.2, k = 0.2, E = 4.4673, L = 0.1. Because of the large 
charge of the test particle the orbital cone has a small opening angle. 



(a) MBO. (b) Projection onto the x-z-p\anc The two cir- 

cles denote the horizons r+ and f — ■ The po- 
tential barrier at r < r~ prevents the particle 
from falling into the singularity. 




(c) BO. (d) Projection onto the s-z-planc. The two cir- 

cles denote the horizons f+. and f _ . 



Figure 8: MBO and BO with parameters Q = 0.45, G = 0.1, q = 0.0, g = 4.0, k = 5, E = 0.9167, L = 0.01. The sphere in the 3d 
plots shows the horizon r_|_. 



(a) TEO. The sphere shows the horizon f+ . (b) Projection onto the a:-z-plane. The two cir- 

cles denote the horizons r+ and f _ . 



Figure 9: TEO with parameters Q = 0.45, G = 0.1, q = 0.0, g = 4.0, k = 5.0, E = 1.0001, L = 0.01. The particle crosses both 
horizons, is reflected at the potential barrier and emerges into a second universe, where it moves out to infinity. 



4.6 The observables 



The Reissner-Nordstrom space-time is believed to be of little astrophysical relevance, because a charged 
black hole would lose its charge quickly by interactions with matter around it. Nevertheless, the study 
of charged black holes can help to understand more general space-times with sources immersed into 
spherically symmetric matter backgrounds. 

For test particle motion one can define frame independent observables such as the perihelion shift for 
bound orbits, the light deflection for escape orbit, the deflecton angle for flyby orbits, or the Lense-Thirring 
effect. We follow the lines of [12] to calculate possible observables. 

We consider here bound orbits of type E and D from the Table [TJ The f-motion is periodic with a 
period 

dy 



UJf 



dr 



(40) 



where e± and ei are the zeros of Ps{y) related to r m - n 

UJr ' 

The period of the ^-motion is given by 

w# = 2 / -= = -2 



and 



The corresponding orbital frequency is 



and the corresponding frequency by 



On 



271 



2tt 



(41) 



The secular accumulation rates of the angle ip and the time t are given by: 



Y v = — 



2 



i-e 



(42) 



and 



T = —IE 
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UJf 
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UJf 



"= 7= + A a / 



A r y/R 



7=2 

■+ 
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(43) 



where I~ and I~ defined in the equations ((37]) and ((38]) are evaluated at j ei corresponding to the root e 
i = 1,2. The orbital frequences and O v are then given by: 

2n 1 _ 2vr 1 

CJf 1 W,J 1 1 



(44) 



The perihelion shift and the Lense-Thirring effects are defined as differences between these orbital fre- 
quences 



^perihelion 
AL cnsc _Thirring 

As expected, there is no Lense-Thirring effect. 

5 Conclusions and Outlook 

In this paper we presented the analytic solution of the geodesic equations of charged test particles in the 
Reissner- Nordstrom space-time in terms of the Weierstrafi p, a and £ functions. The derived orbits depend 
on the particle's energy, angular momentum, electric and magnetic charges, Carter constant and on the 
parameters of the gravitating and charged source. We discussed the general structure of the orbits and 
gave a complete classification of their types. The orbits of charged test particles lie on a cone similar to the 
motion in Taub-NUT space-times with gravitomagnetic charge. This happens because of the non-vaning 
A g cost? term in the equations ([8j and ©. Thus, the motion of only electrically charged particles around 
an only electrically charged source or the motion of only magnetically charged particles around an only 
magnetically charged source reduces to motion in a plane. However, this is not true for the motion of test 
particles with both types of charge around a charged source or the motion of charged particles around a 
source with both types of charges (as in the examples considered in this paper). 

We showed that for large charge of test particles bound orbits behind the inner horizon are possible. 
So, a test particle moves in a space-time, that has horizons which hide the singularity, but for the motion 
of the particle the singularity is not hidden. Such an effect is not present for neutral particles. 

Finally we remark, that as compared to the geodesies in the field of a gravitomagnetic monopole |12j , 
the geodesies in the Reissner-Norstrom space-time are complete, because the space-time can be the ana- 
lytically continued (without destroying some of its essential properties [16]). Thus at least theoretically, a 
particle crossing the horizons may emerge afterwards into another universe, following a many-world orbit. 



- ? 



(45) 
(46) 
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A Integration of elliptic integrals of the third kind 

f v dv 

We consider an integral of the type I\ = I — — s . This integral is of the third kind because the 

Jv in P( v ) ~ Pi 

function fi(v) = (p(u) — pi) -1 has two simple poles v\ and i>2 in a fundamental parallelogram with 
vertices 0, 2wi, 2wi + 2u;2, 2w2, where 2lo\ and 2a_>2 are fundamental periods of p(v) and p'(v). 
Consider the Laurent series for the function /i around vi 

fi{v) — a_i i j(w — «i) 1 + holomorphic part , (47) 

and the Taylor series of /j -1 about Vi 

fr\v) =p' («,-)(« - Vi) + 0(v 2 ) . (48) 

Comparing the coefficients in the equality 1 = fi{v)fi(v) where 

1 = fi(v) (p'(vi)(v - Vi) + 0(v 2 )) = o_i, iP '(wi) + 0(v 2 ) (49) 



yields a_i, = —, — r. Thus, the function fi(v) has a residue ,} ; in Vi 



The Weierstrafi £(u) function is an elliptic function with a simple pole in and residue 1. Then the 
function A\ = fi(v) — X)i=i ^p'^v-) 1 ^ s an elliptic function without poles and therefore a constant [26] . 
which can be determined from /i(0) = 0. Thus, 

h(v)=± C{V - V ^ ] , (50) 
j~[ P{ v i) 

here p'{v2) — p'i^ujj — v{) — — p'(vi). Applying now the definition of the Weierstrafi cr-function 
J v C( v )dv = loger(u) — \oga(vi n ) upon the integral 1\ we get the solution 

U = £ AW* = t ^ (««,)(. - .,») + '"s^^) - (51) 



B Integration of elliptic integrals of the type I2 = f v . 



dv 



in (p(«)-_p 3 ) 2 

We consider the Laurent series of /2(f) and the Taylor series of f^iv) around v t for /2(f) = 

(pO) - Ps) 

f2(v) — a^2,i{v — Vi)~ 2 + a_i^(w — i>i) -1 + holomorphic part , (52) 
fi\v) = ( y p'(v l )(v~v l ) + lp / '(v l )(v-v l f + 0(v 3 ^ 

= (p'K)) 2 (« - «0 2 + P'(«i)p>i)(« - «0 3 + 0(« 4 ) • (53) 

The function /2(f) has poles of second order in v\ and i>2 such that /2(wi) = P3 = j '2(^2) • Comparison of 
the coefficients in 

1 = / 2 («)((p>i)) 2 («-^) 2 + P>i)p>i)(«-^) 3 +0(« 4 )) 

= a_2,i (p'(^)) 2 + [v - Vi) a-i^(p'(vi)f + a^2,ip'{vi)p"(vi) +0(v 2 ) (54) 

yields 

a_2 i = 9 , a-i,i = o . (55) 

(p>0) 2 (p>0) 3 

The function p(v) possesses a pole of second order in v = with residuum and the Laurent series of 
p begins with v~ 2 . Then the Laurent series of a_2,ip(i> — Vi) around Vi begins with a—2,i(v — Vi)~ 2 which is 
similar to the first term in (|52[) . The Laurent series of a— i,iC(u — u<) around Vi begins with 0-i,i(u — Vi) . 

Thus, the function A 2 = /2(f) — ( — — ~ — ll J has no poles and is constant [55] and 

i=i \ (P'K)) (P>i)) / 



can be calculated from /2(0) =0: A2 = — I — - 



p{vi) p"(vi)C(vi) 



^V(p>*)) 2 (p>0) 3 



form: 



Using of f v p(v)dv = —((v) + ((v in ) and the definition of the cr-function the integral I 2 takes the 



r v 2 
j f2(v)dv = A 2 (v - v^) - ^2 

J V\ n ■ 1 



C{V ~ Vi) - C(vin ~ Vi) + — — — log — 

p'(yi) (?{Vin-Vi) 



(p'K)) 2 



(56) 
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